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Abstract 

We study the combined effects of vacuum polarization, relativity, 
Bremsstrahlung, and atomic polarization in nuclear reactions of astrophys¬ 
ical interest. It is shown that these effects do not solve the longstanding 
differences between the experimental data of astrophysical nuclear reactions 
at very low energies and the theoretical calculations which aim to include 
electron screening. 
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I. INTRODUCTION 


Understanding the dynamics of fusion reactions at very low energies is essential to un¬ 
derstand the nature of stellar nucleosynthesis. These reactions are measured at laboratory 
energies and are then extrapolated to thermal energies. This extrapolation is usually done 
by introducing the astrophysical S-factor: 

a {E) = ^ S{E) exp [-27rp(E)] , (1.1) 

where the Sommerfeld parameter, r]{E), is given by r]{E) = ZiZ 2 e^/hv. Here Zi^ Z^-, and v 
are the electric charges and the relative velocity of the target and projectile combination, 
respectively. 

The term exp(—27r?7) is introduced to separate the exponential fall-off of the cross-section 
due to the Coulomb interaction from the contributions of the nuclear force. The latter is 
represented by the astrophysical S-factor which is expected to have a very weak energy de¬ 
pendence. The form given in Eq. (2.1) assumes that the electric charges on nuclei are “bare”. 
However, neither at very low laboratory energies, nor in stellar environments this is the case. 
In stars the bare Coulomb interaction between the nuclei is screened by the electrons in the 
plasma surrounding them. A simple analytic treatment of plasma screening was originally 
given by Salpeter |]l|. In most cases of astrophysical interest Salpeter’s treatment still re¬ 
mains to be a sufficient approximation In the very low energy laboratory experiments 
the bound electrons in the projectile or the target may also screen the Coulomb potential 
as the outer turning point gets very large (> 500 fm). As experimental techniques improve 
one can measure the cross section in increasingly lower energies where the screened Coulomb 
potential can be signihcantly less than the bare one. This deviation from the bare Coulomb 
potential should manifest itself as an increase in the astrophysical S-factor extracted at the 
lowest energies. This enhancement was indeed experimentally observed for a large number 
of systems [^|-§]. The screening effects of the atomic electrons can be calculated [^] in the 
adiabatic approximation at the lowest energies and in the sudden approximation at higher 
energies with a smooth transition in between 

In the adiabatic approximation one assumes that the velocities of the electrons in the 
target are much larger than the relative motion between the projectile and the target nuclei. 
In this case, the electronic cloud adjusts to the ground-state of a “molecule” consisting 
of two nuclei separated by a time-dependent distance i?(t), at each time instant t. Since 
the closest approach distance between the nuclei is much smaller than typical atomic cloud 
sizes, the binding energy of the electrons will be given by the ground-state energy of the 
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Zp + Zt atom, B[Zp + Zp). Energy conservation implies that the relative energy between 
the nuclei increases by Ue = B{Zp) — B{Zp + Zp). This energy increment increases the 
fusion probability. In other words, the fusion cross section measured at laboratory energy 
E represents in fact a fusion cross section at energy E + Ue, with Ue being called by the 
screening potential. Using eq. (1.1), one gets 


a {E -\- Ue) = exp 




a(E), 


( 1 . 2 ) 


where one assumes that the factor S{E)/E varies much slower with E, as compared to the 
energy dependence of exp [—27ir]{E)] . 

The exponential factor on the right-hand-side of eq. (1.2) is the enhancement factor 
due to the screening by the atomic electrons in the target. For light systems the velocity 
of the atomic electrons is comparable to the relative motion between the nuclei. Thus, a 
dynamical calculation has to be done for the effect of atomic screening [Q. However, the 
screening potential Ue obtained from a dynamical calculation cannot exceed that obtained in 
the adiabatic approximation because the dynamical calculation includes atomic excitations 
which reduce the energy transferred from the electronic binding to the relative motion. 

Contributions from the nuclear recoil caused by the atomic electrons are expected to 
further increase the screening effect for asymmetric systems P,ITn 


In almost all the cases 


observed screening effects are found to be equal to or more than the theoretical predictions. 
Recently including improved energy loss data for atomic targets is shown to lead agreement 
between theory and data 110 , however the situation is still not resolved for molecular and 
solid targets. Electron screening enhancement was not observed for the heavier symmetric 
system ^Hei^He, 2p)'^He which is expected to have about 20% enhancement at the 
energies studied. Recent measurements have not yet clarihed the effects of electron 
screening in this reaction. A mechanism which reduces the screening enhancement for this 
system (and possibly for other systems with large values of ZxZ 2 and the reduced mass) 
seems to be needed. 

In this article we show that the contributions from the polarization of the vacuum, 
relativity, Bremsstrahlung, and atomic polarization cannot achieve this task. The motivation 
for this investigation is that Ue/E, appears multiplied by a large number, r]{E), in the 
exponential factor of eq. (1.2). One needs only a small value of Ue to obtain a sizeable 
enhancement factor; typically Ue/E ~ 0.001 . The effects of the vacuum polarizability 
was previously investigated in Ref. |]^ for elastic scattering below the Coulomb barrier 
and in Ref. for subbarrier fusion reactions using the formalism developed by Uehling 
Effects of vacuum polarization in C scattering at 4 MeV was subsequently 
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experimentally observed [Q. Other small effects in elastic scattering at low energies have 
also been studied by several authors |T9|-|^. They have also been studied in the context of 
astrophysical reactions in refs. However, to our knowledge, the other effects have 

not been studied. 

In section 2 we study the effects of vacuum polarization, relativity, Bremsstrahlung, and 
atomic polarization for the astrophysical reactions listed in ref. |^, and for which a set of 
“experimental ” values of screening energies Af/g are given. These experimental values are 
chosen so that the eq. (1.1) reproduces the enhancement of the fusion cross sections at very 
low energies. In section 3 we present our conclusions. 


II. SMALL EFFECTS IN THERMONUCLEAR REACTIONS 


To calculate the fusion cross section corrections we use for simplicity the (s-wave) WKB 
penetrability factor 


P{E) = exp 




( 2 . 1 ) 


where p{r) is the (imaginary) particle momentum inside the repulsive barrier. The corrected 
fusion cross section is given by 


a= ac.R, (2.2) 

where ac is the pure Coulomb repulsion cross section, and R = Pc+a{E)/P{E), with a = 
{scr, VPol, rel, Brems, At}, are the corrections due to atomic screening, vacuum polarization, 
relativity, Bremsstrahlung, and atomic polarization, respectively. 

The atomic screening effect is calculated using |p(r)| = where E 

is the relative energy between the nuclei. The atomic screening potential, assumed to be a 
constant function of r (valid for r -C = 0.529 A), is given by Ug. 


A. Vacuum Polarization 


Vacuum polarization increases the electromagnetic potential between two like charges. 
Like the Coulomb potential itself, the increase due to vacuum polarization is also propor¬ 
tional to the product of the charges [|^ . Vacuum polarization contribution increases almost 
exponentially as the two charges get closer. The Coulomb interaction is smaller for asym¬ 
metric systems than for symmetric systems of comparable size. On the other hand, the 
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nuclear force tends to extend farther out for asymmetric systems because of the extra neu¬ 
trons. Consequently for asymmetric systems the very tail of the nuclear force can turn the 
relatively weak Coulomb potential around to form a barrier at a considerable distance from 
the nuclear touching radius. For symmetric systems, however, the location of the barrier is 
further inside where the vacuum polarization contribution is stronger. We show that the 
resulting increase in vacuum polarization is nevertheless not sufficiently large to make an ap¬ 
preciable contribution to the extracted astrophysical S-factor. For light symmetric systems 
with small values of this effect should be negligible. Indeed, for the pp reaction the 
vacuum polarization contribution was shown to be very small Similarly the measured 
S-factor for the d{d,p)^H reaction [§] agrees well with theoretical calculations of atomic 


screening [21|. On the other hand one may expect that already for the ^He{^He,2pYHe 


reaction the increase in the potential due to the vacuum polarization could be large enough 
to counter the decrease due to electron screening. We show that this is not the case. 

The vacuum polarization potential is according to Uehling given by 


^ 1 ^ 26 ^ 2a (2r\ 

Vpjr) = —— - / - 


(2.3) 


where a = 1/137 is the hne structure constant, and Ae = 386 fm is the Compton wavelength 
of the electron. The function I (x) is given by 


I{x) = / 


oo 1 ,/f2 — 1 


As shown by Pauli and Rose this integral can be rewritten as 


poo 

I{x) = a{x)KQ{x) + /3{x)Ki{x) + 'y{x) / Ko{t) dt , 

Jx 


(2.4) 


(2.5) 


where 


7(x) = ^x{l + 


I3{x) = -^a:(l + 
with X = 2 r/Ae . 


( 2 . 6 ) 


In Ref. [T^ it was shown that the modihed Bessel functions Kq and Ki as well the integral 
over Kq can be expanded in a very useful series in Chebyshev polynomials which converge 
rapidly and for practical purposes only a few terms (?« 5 — 10) is needed, allowing a very 
fast and accurate computation of the Uehling potential. 

In Figure 1 we plot the Coulomb potential and the vacuum polarization potential for 
Z 1 Z 2 = 1. Both the Coulomb potential and the screening potential scale with the product 
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Z\Z 2 - However, the vacuum polarization potential has a stronger dependence on the nuclear 
separation distance. 

The limits of the integral are the nuclear radius, where the nuclear fusion reaction 
occurs, and the classical turning point in the Coulomb potential, Rc = Z\Z 2 ^ jE'^ where 
E' = E + Ue- At very low energies the inferior limit is not important when vacuum 
polarization is neglected (the exponential factor in Eq. (1.2) can be obtained with Vpoi = 0, 
and Rn —>■ 0, in Eq. (2.1)). However, since the vacuum polarization potential has a strong 
dependence on the nuclear separation distance, being much stronger at shorter distances, 
its effect is very much dependent on the choice of this parameter. For all reactions with the 
deuteron we use the “deuteron radius”, Rn = 4.3 fm, corresponding to an average distance 
value associated with matrix elements involving the deuteron. For the other reactions we 
use Rn values given in the third column of Table I. In the 4th row we show the ratio 
between the penetrability factor through the Coulomb barrier and the penetrability factor 
including atomic screening, Pc+scr{,E)/Pc{E). In the 5th row we show the effect of vacuum 
polarization, Pc+vPoi{E)/Pc{E). The energy E chosen is the lowest experimental energy 
for each reaction. The atomic screening corrections Ug were calculated in the adiabatic 
approximation, given by the differences in electron binding energies between the separated 
atoms and the compound atom ||^. We see that the effect of vacuum polarization is small, 
but non-negligible for some reactions. Moreover, it increases the discrepancy between the 
value of the screening potential required to explain the experimental data and the theoretical 
calculations of this potential as illustrated in Table I. 


B. Relativistic effects 


A classical Hamiltonian may be written which contains relativistic effects to hrst order 
in 1/mc^. This Darwin Hamiltonian takes the following form in the center of mass system 


E = 


p 


2mc 


+ 


ZpZre^ 



(2.7) 


where mo = mprriT /(mp + mp) is the reduced mass, and pr is the radial component of the 
relative momentum. In a head-on collision, 

2mo r 8c^ \mp rrip) mprriTC^ r 

The solution of this equation yields, for Rn <r < Rc, 
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\p\ = m 


- 1/2 


(a + 7) + \J{a + 7)^ + 4/9 {Vc - E) 


1/2 


a = 


l/2mo , /9 = (^1/mp + 1/m^ /8c^ , 7 = ZpZxe^/ (mpfriTC^r^ 


(2.9) 


The correction due to relativity is given in the 6th column of Table I. Although the 
correction in the momentum |p| is of the order of 10“®, the penetrability is enhanced by an 
amount of order of 10“^ as compared to the penetrability with only the Coulomb interaction. 

C. Bremsstrahlung 

The energy emitted by Bremsstrahlung per frequency interval du and solid angle element 
dQ is 

|2 


U 


dEpri^) = dudVt— 






( 2 . 10 ) 


where k = uj/c and the direction of observation n = k/Zc. The sum goes over the charges 
positions r^, and velocity \j of the moving particles. In the long-wavelength approximation 




( 2 . 11 ) 


and 


dEpriuj) = dudfl Ap 


UJ 


d(a;) X n —^Q(a;) x h -|- 


( 2 . 12 ) 


where, in the center of mass system, with relative positions and velocities given by r{t) and 
v(t), respectively, and 

d(t) = fie\{t) , 

Q{t) = f 2 e^ [n.r(t)] v(t) , 
f _ 4A-1 ( 

[a^ ’ 

ApAp 


In a head-on collision 


Ap — 


Ap -\- Ap 


( 2 . 13 ) 


r = ao(cosh,^ -|- 1) , 
t = a;o(sinh,^ -h 0 > 


Oq — 


ZpZpe^ 
2E 


Qjq 

Cdo = — • 
V 


( 2 . 14 ) 
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d(ci;), Q(ci;) are the Fourier transforms of d(t) and Q(t), respectively, which can be calculated 
analytically. The hnal result, after an integration over hi, is 


dEsriuj) 4 /V 

—--= -hujoa - 

duj TT Vc 


f;lh + fl A% ft; 


(2.15) 


where 


fci = f [KU^f , ft2 = ^ e-” , 

2 ; = u/uo , 

K..(y) = lf t;os( 4 ) d( , K^M = . (2.16) 

The above results give the Bremsstrahlung due to the incoming branch of the trajectory 
only, since we are interested in the energy loss until the fusion occurs. The result for the 
full trajectory, including the outgoing branch is obtained by replacing the lower limit of the 
integral in eq. (2.16) by —oo. 

The frequencies 2 : -C 1 (cu cuo) dominate the spectrum and we can replace the functions 
Kiz{z) and K'-^{z) by their approximate values for low ; 2 . The integration of eq. (2.15) over 
u is straightforward, and we get 

Ali (“) 

where ({2, 2) = 0.64493... is the Riemann’s Zeta function for a particular value of its argu¬ 
ment. 

The results of energy loss by Bremsstrahlung are given in the 7th row of table 1, where 
Rbv = Pc+Br/Pc is calculated by using eq. (1.2), with Ue replaced by {—Esr)- It is larger for 
the systems with a large effective dipole charge /i, since the quadrupole radiation is smaller 
by a factor (v/c)^. However, even for the systems with a larger /i, the Bremsstrahlung 
correction is of order of 10“^. 


D. Atomic polarizability 

The virtual excitations of the atomic electrons in the target yields in second order per¬ 
turbation theory a polarization potential given by 


u. = - E 


|(0|U(r.R)|n>|'- 

En — Eq 


(2.18) 


where 



Vc{r, R) = E if > r} approx.) , (2.19) 

ZpeVio(r)|^^/f (dipole approx.) , (2.20) 

The first equation is valid in the monopole approximation and the second equation is valid in 
the dipole approximation, in a head-on collision, r* are the positions of the atomic electrons, 
and R is the distance between the atomic nuclei. 

Using hydrogenic wavefunctions and considering only the atomic polarization arising 
from the transitions from the ground state, ^nim = "hioo, and the s-state, <hnZm = *^ 200 , we 
get 

(Oioo V”"‘(r.R)l<l>2oo> = , (2,21) 

27 Qq 

where y = ?>ZTR/‘2ao and 

/(x) = (1 + X) exp(-x) . (2.22) 



For the dipole excitations, considering transitions from the ground state, ^nim = ‘hioo, and 
the p-state, ^nim = *^ 210 , we get 



8 Zj'Zp 

9a/2 Oo 


^(x), 


(2.23) 


where 


^(x) 


8 - (S + 8x + 4x^ + x^) exp(-x) 


(2.24) 


For an estimate of the effect of atomic polarizability, we will assume that the contribution 
of virtual excitations of the s- and p-orbit are the most relevant. For a hydrogen atom 
E 2 — El = 10.2 eV. Using this value in eqs. (2.21-2.24), also for heavier atoms, we find 
the corrections due to atomic polarizability presented in the two last rows of Table I. We 
see that atomic polarizability is more important for monopole excitations. The reasoning 
here is the same as in the use of the adiabatic approximation for the effect of electron 
screening: the monopole held of the combined {Zp + Zp) atom is stronger than the monopole 
held of the Zp atom. Thus, the contribution of the monopole term dominates over other 
multipolarities. However, its ehect on the fusion cross section is still small. This agrees with 
the hypotheses used in the dynamical calculations [^j of atomic screening that ehects due to 
atomic excitations, and particularly for virtual excitations, are small and can be neglected. 


9 









E. Conclusions 


In conclusion, we have shown that the vacuum polarization, relativistic corrections, 
Bremsstrahlung, and atomic polarization contributions to the astrophysical S-factor never 
exceed a few percent, but may be signihcant in extrapolating the measured S-factor to 
lower energies. Although these contributions are not comparable to that of sub-threshold 
resonances and electron screening, they represent some of the many factors that may con¬ 
tribute to the weak energy dependence of the S-factor. Nuclear polarization effects, were 
not included, since they are much smaller than effects due to atomic polarization, for light 
targets. 

Vacuum polarization effects are the most important from all small contributions, and 
sensitive to the inner turning point of the potential barrier, hence to the diffuseness of 
the nuclear potential employed. Although the energy needed to create a virtual e’''e“ - 
pair is much larger than atomic excitation energies, the magnitude of its effect is greatly 
compensated by its large matrix element (due to the large overlap of the electron and 
positron wavefunctions), contrary to the atomic polarization cases. For the same reason, 
nuclear polarization and excitation should be neglected. 

We have shown in this work that none of the corrections beyond the effect of atomic 
screening can explain the missing enhancement of the fusion cross sections in atomic target 
experiments. As suggested in ref. [^, this effect might well be due to a wrong assumption 
on the dependence of the stopping power on the beam energy. 
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Reaction 

-S'min 

[keV] 

bv] 

Rn 

[fm] 

1 Rscr 

1 — Rvpoi 
[x 10-2] 

1 Rrel 
[x 10-3] 

1 — Rpr 
[X 10-3] 

1 - Rir 

|x 10-=] 

_ j^mon 

[x 10-3] 

Did,p)T 

1.62 

20 

4.3 

0.164 

-0.95 

0.17 

0.54 

1.01 

0.246 

^He{d,pYHe 

5.88 

119 

4.3 

0.331 

-1.60 

0.47 

1.12 

0.39 

0.314 

D{^He,pYHe 

5.38 

113 

4.3 

0.364 

-1.58 

0.47 

1.00 

2.11 

0.357 

^HeYHe, 2pYHe 

25 

292 

3.0 

0.196 

-3.14 

1.75 

0.58 

0.35 

0.321 

^Li{p, aYHe 

10.74 

186 

3.0 

0.258 

-1.82 

1.07 

1.36 

0.30 

0.360 

'^Li{p,aYHe 

12.70 

186 

4.3 

0.198 

-1.88 

1.04 

1.28 

0.17 

0.284 

^Li{d, O-YHe 

14.31 

186 

3.0 

0.218 

-2.32 

0.72 

0.71 

0.15 

0.313 

HYli.afHe 

10.94 

186 

3.0 

0.250 

-1.82 

1.07 

1.23 

2.55 

0.350 

HCLi,aYHe 

12.97 

186 

4.3 

0.191 

-1.88 

1.04 

1.17 

1.42 

0.275 

DYLi,aYHe 

15.89 

186 

3.3 

0.184 

-2.34 

0.71 

0.35 

0.91 

0.262 

^^B{p,aYBe 

18.70 

346 

3.3 

0.376 

-2.38 

2.03 

1.45 

0.57 

0.758 

^^B{p,aYBe 

16.70 

346 

2.0 

0.462 

-2.36 

2.00 

2.13 

0.85 

0.906 


Table Caption: 

Lowest experimental energies, -Emin, energy corrections |K| due to the screening by the 
atomic electrons, t/g, nuclear radii, and correction factors for the nuclear reaction: (a) due 
to atomic screening, 1 — Rscr, (b) vacuum polarization, 1 — Rypoi, (c) relativity, 1 — Rreu 
(d) Bremsstrahlung, 1 — Rpr, and (e) atomic polarization, 1 — Rat- 


Figure Caption 


Fig. 1 - Comparison between the Coulomb potential and the vacuum polarization poten¬ 
tial as a function of the nuclear separation distance for = 1. The vacuum polarization 
potential has been multiplied by a factor 1000 in order to be visible in the same plot. 
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